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Abstract. The work on a new fully variational model of average-atom in 
quantum plasmas using a numerical code called VAAQP is reported. A brief 
description of the code is given. Application to aluminium at solid density 
and temperatures between 0.05 and 12 eV is presented. Comparisons to results 
obtained using other approaches are also shown and discussed. The results prove 
the feasibility of the variational model in the warm dense matter regime. Effects 
of the variational treatment can lead in this region to significant differences with 
respect to existing models. 
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1. Introduction 

In this paper, we present the first application in a numerical code (VAAQP) of a new 
fully variational model of quantum plasma at thermodynamical equilibrium. Models 
of screened ions in equilibrium plasmas with quantum all electrons are useful to photo- 
absorption and equation of state calculations. The existing models [U El [3] are not 
fully variational in the sense that they do not respect the virial theorem. Recently 
a fully variational model respecting the virial theorem has been proposed [H [S]. 
The variational character of the model gives a firm basis to the atomic structures 
calculations and provides electronic pressures that are thermodynamically coherent. 

Among the most interesting applications are cases corresponding to the Warm 
Dense Matter (WDM) regime. Special attention is thus paid to matter at high density 
(ie. solid density) and relatively low temperature (around few eV or below). The main 
problem in such regime comes from the long range Friedel oscillations. 

This paper is organized as follows. Sec. II contains a short review of the 
model. Sec. Ill is devoted to a brief description of the VAAQP code. Some results 
obtained with the code for aluminium (Al) in the WDM regime are presented in Sec 
IV. Differences between models requiring the neutrality of the WS sphere and the 
variational model are discussed in this last section. 
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2. Variational theory of average-atom in jellium 

The VAAQP code uses a new model [4j [5] of quantum plasma at thermodynamical 
equilibrium with bound and continuum electrons treated within the same formalism. 
In this model, the equilibrium is determined by 3 parameters: atomic number Z, 
temperature T and ion density rij. All other quantities, including mean ionization Z* , 
are determined from variational equations using for instance the density functional 
theory (DFT). In the Thomas-Fermi (TF) approximation to the electron density the 
present model leads naturally to the Thomas- Fermi ion-in-cell approach [BJ. 

The starting point of the present model is the cluster expansion [7J [HI [H] from 
which the two first orders are retained: 

f(n i ,Z,T) = f + + (1) 

/o is the free energy per unit volume of an uniform electron gaz with an unknown 
density no, 

no = riiZ* (2) 

In the zeroth order, electron charge density is neutralized by an uniform ion 
background. is the average contribution of one ion immersed in that jellium 

to the free energy density. In the first order, the system is no more homogeneous. The 
cluster expansion leads to: 

</>! = UiJ d\ ^fon+iellium (Jf . ^ z j,. ^ _ /q ^ p) j (g) 

The substraction of the zero order term assures the convergence of the above integral. 
All structure variables, including no are denoted by X. 
The free energies per ion are thus defined as follows: 

771 j? I rr-n _ fo (no, T) 

Fq = F (n , rii, T) = (4) 

n l 

F{X-n u Z,T) = F + AF 1 (5) 

m = J d 3 r ^on+jeUium {X; ^ ^ j,. f) _ /q ^ T) | (g) 

All structure variables X are determined from the minimization of the free energy 
F (X; ni, Z, T) with additional conditions: 

• The first order non-central-ion distribution is not uniform but has the form of a 
cavity with a radius R: 

n d(r - R) (7) 

• The neutrality of all charges, taking the cavity into account, is respected: 

Z + J d 3 r{n(r)-n Q 0{r-R)} = (8) 

• The total electron number per unit volume has the following cluster expansion: 

Zn t =n Q + ni J o9r {n(r) — n 6(r — R)} + ... (9) 
Using Eq. [5] in Eq. [5] leads to the identity: R = Rws 
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The minimization leads to the SCF equations plus a new one which allows to 
determine uq or Z*\ 



where V e i(r) is the electrostatic part of the SCF potential [HE]. 

The main difference between iNFERNO-type models [TJ [3] and the present one is 
the way of determining the mean ionization Z*. In the Inferno type models, it stems 
from the neutrality of the WS sphere whereas in the present one it is obtained from 
the variational principle (jlOp . In this model, the WS sphere is in general non- neutral. 

3. The VAAQP code 

The VAAQP code solves the equations of the variational model (see Eqs. H UTU)) . 
First, solutions to the SCF equations in the four-parameter space (Z, T, n,-, Z*) 
are calculated. Then among them is found such that fulfills Eq. [TU] (see FIG.[T|). The 
code has three options to calculate the electron density: semiclassical (Thomas- Fermi) , 
quantum-non-relativistic (Schroedinger) or quantum-relativistic (Dirac) formalism. 

In the case of quantum calculations, continuum states contributions to observables 
are integrated over energy using a resonnance-catching adaptative-mesh-refinement 
method. Electron density is numerically computed for radius from zero up to an 
asymptotic value where it is matched to the asymptotic form given by the linear 
response theory of an homogeneous dense plasma [101 lllj . Solutions are considered 
valid if they are unsensitive to choice of the matching radius and are not mesh- 
dependent. The matching radius depends mostly on the range of the Friedel 
oscillations and becomes larger as temperature is lower (see FIG. [2]). 

Instead of looking for solution to Eq. [TQl the code can be set up to fulfill other 
conditions as for instance the neutrality of the WS sphere. 




(10) 
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4. Model comparison in case of warm dense aluminium 

The code has been applied in the case of Al at typical condition of WDM ie. solid 
density (2.7 g.cm~ 3 ) and temperatures between 0.05 and 12 eV. In this regime 
Friedel oscillations decay on a scale much larger than the WS radius (decay length 
(2.6^) _1 = 55.2 a at 0.05 eV temperature, WS radius Rws — 2.99 a ). In such 
cases, region of numerical calculations should extend over several tens of Rws- It 
is important to recall that the present approach is based on the assumption that 
the SCF potential decays exponentially due to screening which is only valid at finite 
temperatures [TUl [IT] . 

The exchange-correlation term is taken in the local density approximation (LDA) . 
For this exploratory study, in order to have direct comparisons with the Thomas- 
Fermi-Dirac (TFD) model, Dirac exchange term p~2] [13] was used. Iyatomi and 
Ichimaru finite-temperature exchange-correlation term [14] has also been tested. In the 
temperature regime of interest in this study, the differences between results obtained 
using these two approximations appeared to be relatively small. 

FIG. [3] displays the behavior of the mean ionization Z* that was defined in Eq.[2] 
This definition is the only one justified in the framework of [5] and can be applied 
to any model as it corresponds to the asymptotic value uq of the electron density, 
which is related to the chemical potential. It is important to notice that the notion of 
unbound electrons is not related to any well defined observable of a quantum operator. 
In the quantum approaches other definitions are sometimes used. One of those is the 
value of the electron density at the WS radius. This definition is motivated by the fact 
that it is identical to the previous one in the Thomas-Fermi case. Another one is the 
difference between the atomic number Z and the sum of the bound levels occupation 
numbers. These other definitions will not be considered in this paper. 

Results of four approaches are presented on FIG. [31 The first one is the result 
of the variational model with electron density calculated using the non-relativistic 
quantum formalism. The second one is the result of the model with electron density 
calculated using the non-relativistic quantum formalism but in which the Eq. [10] was 




Figure 2. Plot of the electron density versus radius for Al at 0.25 eV temperature 
and 2.7 g.cm~ i matter density with a detailled view of the Friedel oscillations. 
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Figure 3. Mean ionization Z* as denned in Eq. [2] versus temperature in case 
of Al at 2.7 g.cm~ 3 matter density. Solid line corresponds to variationnal model 
with Thomas-Fermi density. Plus symbols denote neutral WS sphere model with 
Thomas-Fermi density. Squares corresponds to variationnal model with non- 
relativistic quantum density. Circles denote neutral WS sphere model with non- 
rclativistic quantum density. 



The condition based on Eq. QT] is used in the Inferno model. However calculations 
reported here with Eq. QT] do not correspond exactly to the Inferno model as in our 
case the potential outside the WS sphere is not set to be zero. Two calculations using 
the TF formalism are also reported, one with Eq. [TU1 another with Eq. QT] According 
to [HH], in the TF case, these two approaches are stricly equivalent which is confirmed 
by FIG. [3] and FIG. HI 

As it is seen on FIG. El results from quantum calculations with Eq. [10] can differ 
significantly from those obtained using Eq. [TTJ Differences are especially pronounced 
at low temperature and tend to vanish as temperature increases. In the case of Al 
at solid density and at temperature below 2.5 eV, the atomic structure is different in 
the two quantum models. The 3s shell is found to exist whithin the variational model 
whereas it is absent when using the NWS model. These differences are mainly related 
to the quantum behavior of the density outside the WS sphere, namely to the Friedel 
oscillations wich are present in this regime. 

Electronic pressure is calculated from the formula: 



/o = /o + fa ° is the free energy per unit volume of an uniform electron gaz including 
exchange-correlation, /q is the perfect gas free energy for the grand-canonical 



replaced by the condition of the neutrality of the WS sphere (NWS): 




(11) 



P = -fa + n (p, + V xc {n ) + V e i(R W s)) 



(12) 
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Figure 4. Electronic pressure P versus temperature in case of Al at 2.7 g.cm~ i 
matter density. Solid line corresponds to variationnal model with Thomas-Fermi 
density. Plus symbols denote neutral WS sphere model with Thomas-Fermi 
density. Squares corresponds to variationnal model with non-relativistic quantum 
density. Circles denote neutral WS sphere model with non-relativistic quantum 
density. 



In the present study f§ c was set to e xc from [13] . V xc is the exchange-correlation part 
of the SCF potential |1 [5]. 



Eq. [T3 is obtained from the variational theory and applies only to the case where 
Eq. [TU] is fulfilled. The quantum NWS model is not variational and the use of the 
above pressure formula is only given for comparison purposes. In the case of models 
that are not fully variational the electronic pressure is often calculated by numerical 
differentiation of the free energy with respect to the ion density. 

As can be seen on FIG. 0] relative differences in calculated pressures from the 
two quantum models can be about 40% in the low temperature region and decrease 
as temperature increases. Pressures obtained using the quantum NWS model at low 
temperatures seem to be relatively close to those obtained from the TFD model. 

FIG. [31 and FIG. H] both display curves that are continuous with respect to 
temperature. In the case of the quantum variational model, it is worth to stress 
that continuity is preserved despite the fact that the 3s shell is disappearing between 
2 and 2.5 eV. This is due to a careful treatment of continuum resonnance in the 



ensemble, namely in atomic units: 




(13) 




(14) 



VAAQP code. 
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5. Conclusion 

In this paper, a fully variational model of the average atom in quantum plasmas 
is applied for the first time in a numerical code called VAAQP. It is shown on 
some chosen examples that the code is able to calculate variational self-consistent 
equilibrium at high matter density and relatively low temperatures. In this regime, 
Friedel oscillations of the electron density and SCF potential present a real challenge 
from the numerical point of view. 

Effects of the variational treatment is studied by comparison to results from a 
model requiring the neutrality of the WS sphere. The variational model can lead to 
pressures, mean ionizations and atomic structures that differ significantly from those 
obtained using existing models especially in the low temperature region. 

The work on the present model is in progress. Among objectives is a coherent 
approach to the equation of state and the radiative properties of the warm dense 
matter. 
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